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A gauge invariant combination of LLe sleptons within the Minimal Supersymmetric Standard 
Model is one of the few infiaton candidates that can naturally explain population of the observable 
sector and creation of matter after inflation. After the end of inflation, the infiaton oscillates 
coherently about the minimum of its potential, which is a point of enhanced gauged symmetry. This 
results in bursts of non-perturbative production of the gauge/gaugino and (s)lepton quanta. The 
subsequent decay of these quanta is very fast and leads to an extremely efficient transfer of the 
infiaton energy to (s)quarks via instant preheating. Around 20% of the infiaton energy density is 
drained during every infiaton oscillation. However, all of the Standard Model degrees of freedom 
(and their supersymmetric partners) do not thermalize immediately, since the large infiaton vacuum 
expectation value breaks the electroweak symmetry. After about 100 oscillations - albeit within 
one Hubble time - the amplitude of infiaton oscillations becomes sufficiently small, and all of the 
degrees of freedom will thermalize. This provides by far the most efficient reheating of the universe 
with the observed degrees of freedom. 



I. INTRODUCTION 

Primordial inflation [j| has many virtues - it explains 
the large scale homogeneity and stretches the initial seed 
perturbations to the observable scales in the cosmic mi- 
crowave background radiation . In spite of the observa- 
tional successes, a prime question remains regarding the 
microphysical origin of inflation - what is the infiaton? 

There are many infiaton candidates, see e.g. [3[, which 
are fully capable of explaining the cosmic microwave 
background radiation, but it is not clear how they 
would create the observed fundamental particles from the 
inflaton-induced vacuum energy. Reheating after infla- 
tion must generate all the Standard Model (SM) degrees 
of freedom required for the success of Big Bang Nucle- 
osynthesis [|| (for a review on reheating, see e.g. @.). 

Inflation can occur within a hidden sector, where the 
infiaton is a Standard Model (SM) gauge singlet, and 
couples to all of the hidden and visible sector degrees of 
freedom. However, the process of transferring the energy 
from the infiaton and all the hidden sectors to the visi- 
ble sector degrees of freedom is largely unknown due to 
the fact that the couplings are mostly unknown beyond 
the strong, weak and electromagnetic interactions. Fur- 
thermore, whether all the hidden degrees of freedom can 
decay into the visible sector or not is a debatable issue. 1 



In order to answer all these questions, one must embed the hid- 
den sector inflation within a known setup where both UV and IR 
aspects of the theory are well understood. Only one such exam- 
ple is known throughout the literature where a closed string 
modulus within Large Volume Compactifications is responsible 



These issues can be resolved if the infiaton belongs to 
the visible sector, and its mass and couplings to the SM 
fields are well known. It has been demonstrated 0-0] 
that inflation can occur within the Minimal Supersym- 
metric Standard Model (MSSM) and its minimal exten- 
sions, with the remarkable property that the infiaton is 
not an arbitrary gauge singlet. Rather, it is a D-flat di- 
rection in the scalar potential consisting of the supersym- 
metric partners of quarks and leptons. 2 . These models 
give rise to a wide range of scalar spectral indices [9| Jl2| - 
including the whole range allowed by WMAP |2| 3 . 

Moreover, since the infiaton is an MSSM flat direc- 
tion, its couplings to matter fields are known. It is there- 
fore possible to track the thermal history of the universe 
from the end of inflation. 4 After the end of inflation, 



for generating the seed perturbations, and all possible couplings 
of the infiaton to hidden and visible sectors are known. In this 
setup, the majority of the infiaton energy density is driven to ex- 
cite the hidden sector particles as compared to the matter fields. 
One requires stringent and very speculative constraints to over- 
come the problem [g|. 

2 For a review on MSSM flat directions, see [l(], [TJ] 

3 Although it was pointed out in original papers that there ex- 
ists tuning of the MSSM parameters at the scale of inflation to 
maintain the flatness of the potential 043- However in subse- 
quent studies it was noted that this tuning can be ameliorated 
substantially not only at high scales, see Refs. but also 
at scales where squark and selpton masses would be measured at 
the LHC, see Ref. EJ- The initial condition problem has been 
addressed in Refs. |l7l ll8(l . 

4 The parameter space permitting successful inflation is compat- 
ible with supersymmetric dark matter flfi , 19] (and may even 
lead to a unified origin of inflation and dark matter |20|). 



the flat direction starts coherent oscillations around its 
minimum, denoted by <j) = 0, which is also a point of en- 
hanced gauge symmetry Q . Whenever the inflaton passes 
through the origin, the eigenmodes of the fields that are 
directly coupled to the inflaton are excited via a non- 
perturbative phenomenon plTI^ . At values of the in- 
flaton away from the minimum, the same modes become 
heavy and therefore it is kinematically unfavorable to ex- 
cite them. However, they can and will decay to particles 
that are coupled to them (but not to the inflaton). This 
is the process known as instant preheating (2f| . 



This scenario has been briefly discussed in its gener- 
ality in [9]. In this paper, we present a detailed inves- 
tigation for the case that the LLe flat direction (which 
consists of sleptons) plays the role of the inflaton. In 
this case the inflaton passage through the origin re- 
sults in non-perturbative production of the electroweak 
gauge/gaugino and (s)lepton quanta. These will sub- 
sequently decay to (s) quarks very quickly. As we will 
show, draining the inflaton energy via instant preheating 
is quite efficient, and nearly 10% of the inflaton energy 
density gets transferred to the relativistic particles at ev- 
ery zero-crossing 5 . 



However, since the inflaton Vacuum Expectation Value 
(VEV) breaks the electroweak symmetry, all of the 
MSSM degrees of freedom will not reach thermal equilib- 
rium immediately. After about 100 oscillations, the am- 
plitude of inflaton oscillations becomes sufficiently small 
and electroweak interactions become efficient. At this 
point full equilibrium is achieved. Due to the hierar- 
chy between the frequency of oscillations ~ O(TcV) and 
the expansion rate ~ O(GeV), the universe thermalizes 
within one Hubble time after the end of inflation. 



II. THE LLe INFLATON 



We concentrate on the LLe flat direction as an infla- 
ton. 6 Here L and e denote the left-handed and right- 
handed lepton superfields, which are doublet and singlet 
under the SU(2)w, respectively, and carry hypercharge 
quantum numbers —1 and +2, respectively, under the 
U(1)y- The D- and F-fiatness conditions require that 
the family indices of the three superfields be different. 
Without loss of generality, we can choose the indices as 
L\Liez- Then the flat direction background can be pa- 
rameterized as 



e3 



> (1) 



Here L, e are scalar components of the corresponding 
superfields and tp is a complex scalar fields. The inflaton 
field <p will be identified by the real part of jp: jd> = ipn. 
The scalar potential along cf> is given by [ 
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4Oto0 there exists an inflection point at 



(3) 



in the potential. Successful inflation can occur within an 
interval \cj> — 0o| ~ (0o/6OMp) in the vicinity of fa. 

For TTiffi ~ O(TeV), as preferred by weak scale super- 
symmetry, we have 4>o ~ O(10 14 ) GeV, and the Hubble 
expansion rate during inflation is H ln f ~ O(GeV). 

After inflation the inflaton starts oscillating about 
<j> = with a frequency m^. The inflaton potential, see 
Eq. J2]), has two points of inflection, 8 



This paper is structured as follows. We begin by pre- 
senting a general discussion of MSSM inflation focusing 
on the LLe inflaton in section 2. In section 3, we derive 
the inflaton couplings to scalar, gauge and fermion fields. 
We discuss reheating after MSSM inflation and the dom- 
inance of instant preheating in transferring the energy 
away from the inflaton oscillations in section 4. In sec- 
tion 5, we discuss the subsequent stage of thermalization 
of the MSSM degrees of freedom, and give an estimate of 
the reheat temperature. Finally, we conclude our paper 
in section 6. Essential expressions for the LLe inflaton 
and the relevant interactions for the scalar, gauge boson 
and fermion fields are given in the Appendix. 



5 To be clear, our discussions concerns creation of particles and an- 
tiparticles associated with matter fields during reheating. Gener- 
ation of matter-antimatter asymmetry happens at a later stage. 



with the respective potential energies 
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(4) 



(5) 



The potential is tachyonic -2m^ < V"(<f>) < between 
the two inflection points, while for <f> < <f>\ we have V" > 
0. 



6 The other possible candidates are udd (where it, d correspond to 
right-handed squarks 7, 9J), and NH U L (where N corresponds to 
the right-handed sneutrino and H u is the MSSM Higgs doublet 
that gives mass to up- type quarks [§[). We will not discuss these 
cases in this paper. 

7 For a detailed dynamics of inflection point inflation, see |13|| . 

8 Because of the potential is symmetric under <f> — > — (j>, there are 
two other points of inflection at — 0o an <i — </>!• 
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An important point to keep in mind is that because of 
the hierarchy ~ 10 3 7?i n f, the inflaton can undergo a 
large number of oscillations within a single Hubble time 
after the end of inflation. As we will see later, this has 
interesting consequences for reheating in this model. 




which includes 10 real degrees of freedom in total. For 
the background given in Eq. (fl}, the inflaton <f> is 



<P\,R + <f4,R + <f5,R 

x/3 



(7) 



The combination (<px,i + f^j + <f5j)/V3 has a mass 
and is irrelevant for inflation (R, I denote the real and 
imaginary parts of a complex scalar field respectively) . 

The inflaton VEV completely breaks the SU(2)w x 
U{1)y symmetry. This results in four massive real scalars 



V>2,R + <P3,R 



<l2.J - <P3J 



X3 



fl.R - fi,R 



V2 



(8) 



whose masses are obtained from the £)-terms (for detailed 
derivation, see Appendix IVIII A) . From Eqs. (|49I51[) we 
find the following mass terms 



FIG. 1: The inflaton potential around the minimum, to- 



gether with the first approximation, 
x = lO 1/8 <?!>/0o. 



The variable 



III. COUPLINGS OF THE INFLATON 

Any time that the inflaton crosses the origin scalar, 
gauge and fermion fields that are coupled to it are pro- 
duced non-perturbatively. Production of particles that 
have couplings of gauge strength to the inflaton dom- 
inates over that of particles with Yukawa couplings to 
the inflaton. After the inflaton passes the origin it rolls 
back to large VEVs, and the produced particles become 
very heavy. They can therefore decay quickly to particles 
that are lighter than them (i.e. particles with no gauge 
couplings to the inflaton). For the LxL 2 e^ inflaton, these 
include all of the (s)quarks, the Higgs and Higgsino parti- 
cles, and L 3 , ex, e 2 (s)leptons. These light fields interact 
among themselves, which may lead to their thermaliza- 
tion. For a proper treatment of these processes, we first 
need to identify different fields and their couplings to 
each other. 



A. Coupling to scalars and their decay widths 

We focus on the scalars that have gauge couplings to 
the inflaton as they play the dominant role in transferring 
the inflaton energy away. The relevant interaction terms 
arise from the _D-term part of the scalar potential. The 
field content of the L\, L 2 , §3 is 



L x = 





, L 2 = 
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P-3 = <f5, 
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Here gw, gy are the SU{2)w and U(l)y gauge couplings 
respectively. 

The remaining four degrees of freedom that are orthog- 
onal to 4> (and its imaginary counterpart) and x' s are 
Goldstone bosons from breakdown of SU(2)w x U(l)y- 
They are eaten by the Higgs mechanism and give rise to 
longitudinal components of the electroweak gauge fields. 
In the unitary gauge, they are completely removed from 
the spectrum. 

The x particles decay to squarks, the Higgs parti- 
cles, and the L3, ei, ei sleptons. The relevant inter- 
action terms between %'s and these fields are given in 
Appendix IVIII Al From Eq. (|52p. and after summing 
over all decay channels, we find 



r — f — r — 

L Xl ~ 1 Xl — L X3 — 
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Note that the decay rate is proportional to the VEV of 
the inflaton, which sets the mass of x fields. 



B. Coupling to gauge fields and their decay widths 

Couplings of the inflaton to the gauge fields are ob- 
tained from the flat direction kinetic terms (see Appendix 
IVIII Bp . From Eqs. (|53I54[> . and after using Eq. 0, we 
find the following mass terms for the electroweak gauge 
fields 



\g Y <fB»B, , 



where 
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Wx - iW 2 
V2 



Wx + iW 2 
V2 : 



(11) 



(12) 



3 



and Wi tfl and are the SU {2)w and U (l)y gauge fields 
respectively. The gauge fields decay to (s)quarks, Higgs 
and Higgsino particles, and L3, e±, &i (s)leptons. The 
relevant interaction terms are given in Eqs. (|551 15€>|) . The 
total decay widths of gauge fields are found to be 



— MV- — ^w 3 — 



(13) 



^vrV^ ' 16tt\/2 
where we have used Eq. (|TT|) for the gauge field masses. 

C. Coupling to fermions and their decay widths 

Couplings of the inflaton to fermions are found 
Eq. dBDD, see Appendix IVlII CI They result in the fol- 
lowing mass terms for fermions 



V6 
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where $2, ^3, ^4 are Dirac spinors defined in 

Eq. (f6"2"j). The ^P's decay to (s)quarks, Higgs and Hig- 
gsaino particles, and -L3, e%, ei (s)leptons. The relevant 
interaction terms are given in Eq. (|64[) . After summing 
over all final states, we find the following decay widths 
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where we have used Eq. (|T4")) for the mass of \l/'s. 

Note that Xi,2,3 scalars, SU(2)w gauge fields 
W^, W3, and ^1,2,3 fermions have the same mass, 
see Eqs. (19I11I14[) . and the same decay width, see 
Eqs. (|10I13I15[) . Similarly, X4 scalar, U(l)y gauge field 
B, and ^4 fermion have the same mass and decay width. 
This is expected from supersymmetry. The inflaton 
VEV results in supersymmetry conserving masses for the 
scalars, gauge fields and fermions. After the breakdown 
of SU{2)w x U(1)y, the bosonic degrees of freedom in- 
clude four real scalars Xi an d four massive gauge fields, 
and their fermionic partners are grouped into four Dirac 
spinors. 



IV. REHEATING AFTER MSSM INFLATION 

We now discuss the reheating stage after MSSM infla- 
tion and various channels for energy transfer from the 
inflaton to SM quarks and leptons and their supersym- 
metric partners. We will focus on the various ways of 
exciting the light and heavy degrees of freedom. 

After inflation, the inflaton starts oscillating about the 
origin with an initial amplitude 0o — 4"3 ■ There are vari- 
ous mechanisms for particle production from an oscillat- 
ing condensate. Here we discuss those that are relevant 
for the MSSM inflaton. 



A. Tachyonic preheating 

As pointed out before, the inflaton potential is tachy- 
onic V"{(f)) < for ((j) /V3) < <j) < (fio- Quantum fluctu- 
ations of the non-zero inflaton modes are amplified due 
to tachyonic instabilities as the field sweeps this inter- 
val [26|. This effect will go away once the amplitude 
of oscillations (j) drops below <fio / v3, which amounts to 
transfer of 10% of the energy density in the oscillating 
condensate. A question arises that to what extent tachy- 
onic preheating is efficient in reducing the energy of the 
inflaton zero mode by this amount. 

Examination of the potential ([2]) shows that -2mJ < 

V"{4>) < within the [(0 O /V3), <t> ] interval. There- 
fore modes with momentum < k 2 < 2m^ can un- 
dergo amplification because of tachyonic instabilities, 
i.e. rik ~ e 4m *** ; where is the time spent during 
the tachyonic phase [26| . The amount of amplification 
needed to increase the energy density in these modes 
from its initial value 0(m^) to lO~ 1 V((/)o) ~ lO~ 2 ro 2 0o 
is O(10 12 ) GeV 4 . This would require at least O(10) os- 
cillations. 

On the other hand, as shown below, instant preheat- 
ing transfers 0(20%) of the inflaton energy to relativistic 
particles in just one oscillation, which is mainly due to 
the hierarchically larger phase space for resonant particle 
production, see Eq. (ITS]) . This implies that instant pre- 
heating reduces ip below po/^/S after O(l) oscillations, 
after which tachyonic preheating will be irrelevant since 
V" > at all times, see Figs. 1 and 2. 

We therefore conclude that tachyonic preheating ends 
before it can become competitive with instant preheat- 
ing in draining the energy density from the inflaton zero 
mode. 



B. Instant preheating 

The fields that are coupled to the inflaton acquire a 
VEV-dependent mass that varies in time due to the in- 
flaton oscillations. For illustration, we first focus on the 
Xi scalar, see Eq. (jH]). 9 The xi quanta are produced ev- 
ery time the inflaton passes through the origin 25] . The 
Fourier eigenmodes of Xi have the corresponding energy 



\Jk 2 + m\ x + 4m 2 i T 2 gi , 



(16) 



with 4>(t) being the instantaneous VEV of the inflaton 
and m Xl the bare (time-independent) mass of the Xi 



9 Similar arguments and calculations will hold for gauge bosons 
and fermions. Single-crossing occupation numbers in instant pre- 
heating are insensitive to the spin of the field. 
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field, and we have written it in terms of the broad reso- 
nance parameter 



<li 



9w4% 
24ml 



> 1. 



(17) 



and the time r = m^t after the inflaton zero-crossing. 
As is well known, a given mode gets excited when the 
adiabaticity condition is violated u>k ~ This happens 



for modes with k < k n 
the origin, where 



each time the inflaton crosses 



i~* 



-^=gw<l>o = 2ml v^i 



(18) 



Here </>o is the velocity of the inflaton at zero crossing. 10 
The growth of the occupation number of mode k can be 
computed exactly for the first zero-crossing, 



Hfe. 



Xi 



exp 



exp 
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7r(fc 2 + to 2 J 



< i . 



(19) 



The total number density of particles thus produced fol- 
lows 



At the time of decay we have </>(idcc,i) — 0o^dcc,i (as- 
suming tdec,i "C n-0 1 ), where idcc.i ~ T^ 1 denotes the 
time between the zero crossing and Xi decay. After using 
Eq. GDJ), we find 



1 > T dc 



( An 



V 3 5w 
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<7i 1/4 » t,,,! . (22) 



For typical values of m^> ~ 100 — 1000 GeV, we have <f> ~ 
0o — 10 14 — 10 15 GeV, at the beginning of oscillations, 
which results in idcc,i ~ 6 x lO^^m^ 1 . This confirms 
the validity of our approximation and ensures that xi' s 
indeed decay promptly after their production. We also 
note that gw4> 2 (t& eCy i) 3> <fio, which ensures that Xi's are 
non-relativistic at the time of decay, and hence using the 
decay rate in the xi rest frame (|10p is valid. 

The energy density in Xi particles soon after zero cross- 
ing is given by 



Pxi ( T ) 
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where 
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This expression corresponds to the asymptotic value and 
assumes there is no perturbative decay of the produced 
X particles. However, immediately after adiabaticity is 
restored, 



x/6 1 _ 

t>t*,l = \{ y => T 1 >r^ 1 = —=q 1 

V gw9o V2 



1/4 



(21) 



Xi particles can and will decay into lighter particles (i.e. 
those particles that have no gauge coupling to the infla- 
ton). As mentioned earlier, in the case of L\Liez infla- 
ton these are the (s)quarks, Higgs(inos), and L3, ex, e2 
(s)lcptons. The transfer of the inflaton energy to these 
light particles is called instant preheating [251 ] . 

As the inflaton is rolling back the potential, the in- 
stantaneous mass that it induces for the produced xi 
quanta increases, and so does their decay width T Xl . 
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and Ki(z) is the modified Bessel function of the second 
kind, satisfying zK\(z) = 1 + 0(z 2 ). 

p Xl actually decays with the decay rate T Xl , 



Pxi ( T ) = Pxi ex P 



r xi dt 
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with Adec,i = ^4i(Tdcc) = (47r 2 /3<7 2 y) ~ 36.6, which can 
be integrated beyond the time of decay, to give 



1i m t 

p Xl ~7.82 x — exp 
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(26) 



with — (f> 2 /2 the inflaton energy at each zero-crossing. 
Thus the fraction that is transferred from the inflaton 
to xi's, and through their prompt decay into relativistic 
squarks at every inflaton zero-crossing, is given by 



10 Neglecting the expansion of the universe, we have rj>o = y 2V (<j>) 
from the conservation of energy, where rfr is the amplitude of the 



inflaton oscillations. 
mjy(t>, see Fig. 2. 



Note that after a few oscillations, 



Pxi 
P<t> 



0.0067 gfo exp 



Xl 



(27) 



The expressions as in Eqs. (I22I26I27I) also hold for X2, X3 

scalars, W + ,W~ ,Wa gauge bosons, and ^i, ^2, ^3 
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fermions. The reason being that all of these fields have 
the same mass and decay width as xi 7 an d single-crossing 
occupation numbers in instant preheating are insensitive 
to the spin of the held. 11 

On the other hand, the Xi scalar, B gauge field, and 
'J 4 fermion have a different mass and decay rate. The 
corresponding expressions for them are 



= 9 Y <i>o 



and 



1 ^> TdccA = ?n<ptdec,4 
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(28) 
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does not modify signihcantly the inflaton oscillations. 
The x variance induced at instant preheating can be com- 
puted exactly, 



<x 2 



k 2 dk rik 



2n 2 Lo k {r) 



fr x dt 



2tt 3 
x exp — 



Ae A 



A 



2A, 



K X (A)-K (A) 



exp 



(34) 



which decays exponentially with time and never consti- 
tutes any danger for the rate of inflaton oscillations, e.g. 
at the time of decay it is less than a few percent of rni. 



which gives the same expression for the number density 
and energy density in Xi particles, but with q = q± and 
Aicc,4 = (8^79^) ~ 24.4, 
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4> 3/2 
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and thus 
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(31) 



(32) 



is significantly larger than that coming from each of the 
Xi,2,3 fields. The expressions in Eqs. (|29I31I32I) also hold 
for the B gauge boson and \&4 fermion, which have the 
same mass and decay width as X4- 



D. Rate of Energy Transfer 

To calculate the rate of energy transfer via instant 
preheating, we must add the contributions from all 32 
degress of freedom: 4 from scalars (xi,2,3,4)j 12 from 
gauge bosons (VK ± , W3, B), and 16 from fermions 
(^1,2,3,4)- The fraction of transferred energy through 
Xi,2,3 scalars, W ± , W3 gauge bosons, and ^1,2,3 fermions 
(24 degrees of freedom in total) follows Eq. (|27|). while 
that through the %4 scalar, B gauge boson, and ^4 
fermion (8 degrees of freedom in total) is given by 
Eq. (|32|) . Accordingly, the total energy transferred to 
(s)quarks via instant preheating will be p re i = 24p xi + 
8/? X4 . For gy ~ gw ~ 0.6, this results in 



Prel 



10.6% (per zero — crossing). 



(35) 



Note that this fraction is independent of the amplitude of 
oscillations. Therefore the energy density of the inflaton 
and quarks/squarks after N oscillations of the inflaton 
field respectively are 



= 0.79^0 : 



p m} = (1 - 0.79 w )p 



(36) 



C. Backreaction of \ particles on the inflaton 



where po ~ <I>q/2 is the inflaton energy density at the 
beginning of oscillations. 



As we will show here, the effective inflaton mass in- 
duced by the produced x quanta, 



= m l + q9wY1 &i ) + 2 9 ^ ^> > ( 33 ) 

i=l 



11 The fraction of the inflaton energy transferred via Yukawa cou- 
plings is given by Eq. II27I I. with the gauge coupling gw replaced 
by a Yukawa coupling. It is seen from the parametric form of 
this expression that the fraction becomes smaller for smaller cou- 
plings. This is why we can safely ignore all fields that have only 
Yukawa couplings to the inflaton. 



V. THERM ALIZATION 

In this section we discuss thermalization of particles 
produced via instant preheating and establishment of 
thermal equilibrium among all of the MSSM degrees of 
freedom. This is when reheating completes and one can 
assign a reheating temperature to the universe. 

A. Thermalization of (s)quarks 

As we discussed, each time that the inflaton crosses 
the origin, about 10% of its energy density is transferred 
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FIG. 2: The inflaton oscillations around the minimum of the 
potential after inflation. The red line corresponds to the infla- 
ton evolution without including the particle production. The 
blue one takes into account the non-perturbative production 
of the x quanta, and their perturbative decay to relativistic 
particles, which drain energy from the inflaton. 



into relativistic (s)quarks. The produced (s)quarks do 
not have a thermal distribution, and the question is how 
long it takes for them to thermalize. 

In thermal equilibrium, the number density and en- 
ergy density of relativistic species are given by n t hr — 
(C(3)/7r 2 ) ff>t T 3 and p thr = (tt 2 /3(%*T 4 , respectively, 
where is the number of relativistic degrees of freedom 
and T is the temperature. For a thermal distribution of 



(s)quarks, where = 135, this results in p\(^ 



1/3 



On the other hand, the number density and energy 
density of (s) quarks after the first inflaton zero-crossing 
are given by n ro i = 24n Xl +8n X4 and p le i — 24p Xl + 



X4 ' 



where n Xl , p Xl , n X4 , p Xi are given by Eqs. (|20I27I30I32[) 



1/4 



1.6n^g[ 3 . This im- 



respectively. We then find that p 

plies that the number density of (s)quarks must increase 
by a mild factor in order for a thermal distribution to be 
established. 

The main processes that increase the number density 
of particles are 2 — ► 3 scatterings of (s) quarks (with gluon 
exchange in the i-channel) that result in emission of glu- 

ons of energy ~ pV. 4 (28|. These processes are effective 
because the SU(3)c symmetry is not broken by the in- 
flaton VEV, and hence the gluons remain massless at 
all times. Upon production, each gluon participates in 
the subsequent scatterings, which implies an increasingly 
faster rate for the process. It is therefore suggested that 
the rate for the 2—^3 scatterings can be considered as 
the rate for thermalization of colored fields. This is esti- 
mated to be [H 



thr 



ac 



1/2 ' 
Prel 



(37) 



where ac is the SU(3)c fine structure constant. 

4 4/3 



0.1,00, sec Eq. ([35]). we find that T^ hl . > m f 12 This 
implies that a thermal bath consisting of the colored 
particles indeed forms much earlier than the next zero- 
crossing of the inflaton. It initially carries ~ 10% of 
the inflaton energy density, which results in a tempera- 
ture T ~ O.2(rrt0</)o) 1 / 2 . The energy density, hence tem- 
perature, of the thermal bath grows after each inflaton 
zero-crossing as more (s)quarks are produced via instant 
preheating. 

Because of the large Yukawa coupling of the top quark, 
H u and its Higgsino partner H u are also brought into 
thermal equilibrium with the colored particles. However, 
the thermal bath does not contain (s)leptons because 
they interact with (s)quarks only via SU(2)w x U(l)y 
gauge interactions, which arc suppressed due to the 
breaking of this symmetry by the large inflaton VEV. 13 

Therefore combination of instant preheating and ther- 
malization via SU (3)c gauge interactions provides a very 
efficient mechanism for transferring the energy from the 
inflaton to a thermal bath of colored particles. The trans- 
fer of energy results in a continuous decrease in the ampli- 
tude of oscillations 0, while temperature of the thermal 
bath T keeps increasing. As a result, this thermal bath 
will dominate the energy density after several oscillations. 



B. Thermal masses 



Production of x particles is affected by their mass m x 
according to Eqs. (|20I30I) . For to 2 > the number den- 
sity of x is exponentially suppressed, while m 2 < re- 
sults in an enhancement of n x . Particle production oc- 
curs due to violation of adiabaticity in u)k, see Eq. (|16p . 
each time that the inflaton crosses the origin. A posi- 
tive m 2 > increases u>k and decreases u>k- As a result, 
the adiabaticity condition is violated for a shorter pe- 
riod of time, as compared with the m x — case, which 
suppresses particle production. A negative m 2 has the 



opposite effect and enhances particle production. 

An important contribution to m x is the effective mass 
induced by x coupling to fields that are in thermal equi- 
librium. This results in thermal masses, which can be 
derived from Eqs. (I49I51() : 



VD 9 f (Xi- 



3=1 



2 |d, 



After using the relations p 1G i w (1.6) n,^ and p re i 



12 A more careful treatment that takes various effects (such as mul- 
tiple scatterings) into account results in a rate that is somewhat 
different from that in Eq. (137 6 (we would like to thank Guy 
Moore for clarifying this issue) . This however does not affect our 
results, since the more precise value for r t ji r will still be much 
larger. 

13 The electroweak symmetry is restored over a short period each 
time that the inflaton crosses the origin. However, this time is 
too short for the gauge interactions to bring the (s)leptons into 
equilibrium with the thermal bath. 
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Here we have only kept those fields that are in thermal 
equilibrium, namely the squarks and H u . 

In thermal equilibrium, for any complex scalar $ we 
have ($ 2 ) = T 2 /6 (see, for example [27[). This results in 



m . 



,cff 



'xs.eff 



9 2 yT 2 
12 



rn . 



9 2 yT 2 



12 



(38) 



At first, the tachyonic thermal masses for Xi,2,3 may 
look confusing. As discussed in Ref. [13], the individ- 
ual contributions of scalar fields to the thermal mass of 
another scalar (induced via D-term interactions) can be 
positive or negative. If all scalars are in thermal equi- 
librium, we then have w 2 ff > 0. However, in our case 
only the squarks and H u are in equilibrium. The other 
scalars have electroweak interactions with the thermal 
bath, which are suppressed by the large inflaton VEV. 
The electroweak symmetry is restored for a very short 
period of time when the inflaton crosses the origin. How- 
ever, the time is too short to bring the inflaton and x's 
into equilibrium. 14 

We note that the tachyonic thermal masses of xi,2,3 
are important only at the inflaton zero crossing. As the 
inflaton rolls back to large VEV's, its contribution to u>k 
dominates that from the thermal bath, which prevents 
the latter from inducing instabilities. 

We also note that m x . c g is zero initially (i.e. before the 
first burst of particle production) but increases each time 
that the inflaton crosses the origin and produces squarks 
via instant preheating. We can equate the energy density 
in these relativistic particles with that of the produced x 
particles, 



p Iel = 0.1055, 



30 



.9* 



(39) 



where g* = 172.5 is the number of thermal degrees of free- 
dom (i.e., the colored particles plus H u and H u ). Thus 



m. 



9 2 y ( 3.2 



or 



12 \ir 2 g* 



9Y /3.2 



1/2 



1/2 



HA 12 \9 



1/2 



(40) 



(41) 



Writing these masses in terms of the initial inflaton en- 
ergy we arrive at the expression 



Prcl 



0.058 exp 



\/P4>/Po 



158 



0.048 exp 



\/p<f>/po 



158 



(42) 



for the inflaton energy converted into relativistic species 
every half oscillation, m^Ai ~ 5. 



Coupled evolution of the relative energy 
densities 



We can thus describe the evolution after inflation as 
a coupled system in which the inflaton energy density is 
transferred to the x fields, whose particles promptly de- 
cay into relativistic particles that immediately thermalize 
due to the rapid rate of strong interactions. This leads 
to a set of coupled differential equations, 



r 



(4) 



exp 



\/P4>/Po 



158 



(43) 



Pxi = r 



(i) 



Px4 = r 



(4) 



exp 



exp 



\/p<p/Po 



158 



r xi VP^Pxi ( 44 ) 



x/pIJpo 



158 



r x4\/P^Px4 (45) 



Prcl 



r xi y/P4Pxi + r A 



P<P PX4 



(46) 



where r^ 1 -* ~ 0.0116 m$ denotes the "effective" decay 
rate due to the loss of energy into Xi,2,3 (as well as 
W ± , W 3 and #1,2,3) particles, I^ 4) - 0.0096 m denotes 
the "effective" decay rate due to the loss of energy into 
into X4 (as well as B and #4) particles, see Eq. (|4"2"j) . 
Also, r xi(4) ~ 1.2(1.7) x lO lo m0 denote the "effective" 
decay rate of Xi, 2, 3(4) particles into (s)quarks, in units of 
^Jp~ip oc 4>. The latter decay rate is so large that there is 
simply no time for building up particle occupation num- 
bers (and thus resonant production like in parametric 
resonance) of x fields. 

We have solved the set of differential equations in 
Eqs. ((33143,44,45) for initial conditions p$(Q)/p = 1, 
Pxi(°)/Po = 0, p X4 (0)/p o = and p IC \(0)/p o = 0, and 
plotted the results in Fig. [3] 

We have also performed the integration of the above 
coupled equations assuming that the thermal masses are 
absent (ignoring the exponential terms in Eqs. (42-45)), 
and we see no appreciable difference. 

It is seen from Fig. 3 that 99% of the inflaton energy 
density is transferred within 20 or so oscillations, which 
is in agreement with the expression in Eq. (|36[) . There- 
fore, we conclude that the inflaton energy decays into 
a thermal bath (consisting of colored particles, H u and 
H u ) within O(20) oscillations, long before the expansion 
of the universe has diluted its energy. As a consequence, 
we can assume that reheating is very efficient in this type 
of theories. 



14 Continuous transfer of the inflaton energy to the thermal bath 
via instant preheating lowers the amplitude of its oscillations. 
Eventually, when the amplitude becomes comparable with the 
temperature, electroweak interactions become efficient and all 
degrees of freedom will come to equilibrium. At this point we 
will have rri^ cff > 0, as expected. 



L thr ' 



D. Reheating temperature of the universe 

Even though the colored particles thermalize at a rate 
see Eq. (|3"T)) . which is much faster than m^, all of 
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reheat temperature T r h to the universe at this moment, 
which is given by 



FIG. 3: The evolution of the different energy densities, nor- 
malized to the initial energy density in the inflaton po, as a 
function of time (in units of m^) according to Eqs. (42-45). 
Note that most of the energy of the inflaton, p$, is transferred 
to relativistic particles after a few inflaton oscillations. 



the MSSM degrees of freedom do not reach equilibrium 
as quickly. Full thermal equilibrium is achieved when all 
of the gauge interactions become efficient [29| . As men- 
tioned earlier, the electroweak interactions are initially 
suppressed, because the large VEV of the LLe inflaton 
breaks the SU{2)w x U(l)y symmetry completely. 

The amplitude of inflaton oscillations (f> decreases as 
the energy is continuously transferred via instant pre- 
heating. Eventually, we reach a point when the VEV- 
induced masses of the \ scalars, electroweak gauge 
bosons and \& fermions, see Eqs. (|9lllll4j) . is smaller than 
the average kinetic energy of particles in thermal equilib- 
rium at all times, i.e. 



gw4> 9y4> 



3T. 



(47) 



After using Eq. ([36)l. relation p rcl « p = (tt 2 /30).9*T 4 

"1/2 

(where we take {?* = 172.5), and the fact that <j) , 
we find that it takes N ~ 0(100) oscillations for this to 
occur. 

At this point all of the MSSM particles are kinemati- 
cally available to the thermal bath and all of the gauge 
interactions act efficiently. In particular, the 2 — » 2 scat- 
terings via electroweak gauge interactions happen at a 
rate ~ a 2 EW T, with ctEw denoting the electroweak gauge 
fine structure constant, which is much larger than m^. As 
a result, scatterings of (s)quarks will bring (s)leptons and 
electroweak gauge/gaugino particles into thermal equilib- 
rium, and destroys the residual inflaton condensate, very 
rapidly. 

Therefore, we conclude that all of the MSSM degrees 
of freedom thermalize and reheating completes within 
0(100) oscillations. Note that due to the hierarchy be- 
tween H in { ~ 1O _3 TO0, this happens within a single Hub- 
ble time after the end of inflation. One can then assign a 



T ^=(jy 1/4 ^o /4 ^2xl0 8 GeV, ( 48 ) 

where we have used — 228.75 (all degrees of freedom in 
MSSM) and p = (4/15)m|0g, see Eq. ©. We note that 
this reheating temperature is compatible with the Big 
Bang Nucleosynthesis bounds on gravitino production for 
a gravitino mass m 3 / 2 ~ C(TeV). It is also sufficiently 
high that various mechanisms of baryogenesis may be 
invoked to generate the observed baryon asymmetry of 
the universe. 

We can summarize the relevant time scales in reheating 
after MSSM inflation as follows: 

• Tosc = 27TTOT 1 : the period of inflaton oscilla- 
tions, which separates successive bursts of non- 
perturbative particle production. 

• 70 ~ O(20)t osc : the time scale for efficient transfer 
of the inflaton energy to relativistic particles via 
instant preheating, which can be considered as the 
effective inflaton decay lifetime, 

• Tthr ~ O(100)t osc : the time scale for all observable 
degrees of freedom to come into equilibrium, which 
represents the end of reheating, 

• Tcxp ~ lOOOm^" 1 : the Hubble time right after the 
end of inflation, which sets the expansion rate of 
the universe. 

It is important to note that the time scales involved in 
microphysical processes r osc , r^, r t hr are nicely separated 
and they are all smaller than one Hubble time T eX p- This 
has simplified our analysis of reheating after MSSM in- 
flation. 



VI. CONCLUSION 

In this paper we have studied an inflationary model 
that naturally gives rise to a successful and efficient re- 
heating of the observable sector and excites the MSSM 
degrees of freedom. As a result, it naturally explains the 
creation of matter that makes the structures we see in 
the universe today. 15 We focused on the particular case 
where the LLe flat direction is the inflaton. We identified 
the inflaton couplings to the scalars, gauge bosons and 
fermions and their couplings to other fields that are not 
coupled to the inflaton. We showed that the the former 



15 An alternative scenario is Higgs inflation |3C| . although reheating 
seems to be less efficient in this case [3 ll | ■ due to the large VEV 
of the Higgs-inflaton after inflation. 
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are produced non-perturbatively every time that the in- 
flaton passes through the minimum of its potential. As 
the inflaton rolls back the potential, they become heavy, 
which results in their very fast perturbative decay to the 
latter. 

The combination of these effects results in a very ef- 
ficient transfer of the inflaton energy to relativistic par- 
ticles via instant preheating, which drains about 20% of 
the energy density in a single oscillation. After O(20) 
oscillations virtually all of the inflaton energy density 
is transferred into relativistic particles, and full thermal 
equilibrium in the observable sector will be achieved af- 
ter 0(100) oscillations. This takes less than one Hubble 
time, which provides by far the most efficient reheating 
of the universe with the observed degrees of freedom. 

This study also lays down the foundations for reheating 
in inflationary models where the inflaton couples to all 
the degrees of freedom with the Standard Model gauge 
couplings. Our analysis holds in general at a qualitative 
level, and the quantitative results can be applied with 
proper modifications. 
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VIII. APPENDIX 



Here $i denotes a scalar field with hypercharge yi, and 
T a are the SU(2)w generators. 

Writing L\, L 2 , S3 in terms of the the inflaton </> and 
x's, see Eq. ©, the electroweak Z?-terms read 

D X = ^-^ + Q\T^Q i + Hl d T 1 H u , d + L\T^U 



Do = 



D 3 = ^-^+QlT 3 Q i + Hl d T 3 H a , d + LlT 3 L 3 



<t>X4 xl-xl-xl-xl 

+ \(m 2 - 4N 2 + 2K| 2 ) + (iSil 2 + |e~ 2 |) 2 ) 



\H,, 



(51) 



where i denotes the family index of squarks (we have 
omitted color indices for simplicity). 

The terms that are relevant for the decay of x's to 
squarks, Higgs particles, and L 3 , ei, e 2 sleptons are 



V D 



9w± 
V2 



Xi (Q\T x Qi + HlT x H u + HlT x H d + l|t 1 X; 
X2(Q\T 2 Q l + HlT 2 H u + H\T 2 H d + LIt 2 L 



X3 



[Q\T 3 Q l + HlT 3 H u + HlT 3 H d + L 3 T 3 L ; 



g(|O l | 2 -4|^| 2 + 2|d. i | 2 ) + (|g 1 | 2 + |g 2 |) 2 ) 



-(\H u \ 2 -\H d \ 2 -\L 3 \ 2 ) 



An interaction term fjx<t'*4> between \ an d a massless 
scalar (f> results in a decay rate a 2 /16irM, where a is a 
coupling of dimension mass and M is the mass of %. 



Here we present the essential expressions that are 
needed to derive the inflaton couplings to scalars, gauge 
bosons, and fermions as well as the decay widths of the 
latter. The starting point is the MSSM Lagrangian, 
which can be found, e.g., in [32]. 



Scalar interactions 



B. Gauge field interactions 

Couplings of the inflaton to the gauge fields are ob- 
tained from the flat direction kinetic terms 



Couplings of the scalar fields to the inflaton are gov- 
erned by the electroweak D-terms: 

Vd ^\gw{Dl + D 2 + D 2 ) + \g 2 Y D 2 Y , (49) 



where 



D a = Y,$\T a $ i , Dy = J2lvM$i- (50) 



where 



D^L 2 



a=l 



(9^ + -gyBft - ig w 



a=l 

(0„ - igyBjea ■ (54) 
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Here W 2 ^, and are the gauge fields of 

SU(2)w and U(l)y respectively. 



fermions are given by 



The gauge fields decay to (s)quarks, the Higgses and 
Higgsinos, and the L 3 , ei,e2 (s)leptons. The interaction 
terms that are relevant for the decay to the scalars are 
given by 



C D ^(Ql -W+'^Q^ - d^Q* Uii W+^Q dti ) + h.c. 

+ ^(Ll*W^d,Ll-d,Ll*W^Ll) +h.c. 

+ ^{Hl*W + ^d,Hl-d,Hl*W + ^Hl)+h.c. 

+ ^(H^W+^Hj d^/W+^Hj) + h.c. 

+ ^{Ql^d^ - Q* d ,M^Qd,i) + h.c. 

+ ^{Ll'WgdpLl- LfWgdpLl) + h.c. 

+ f(i/X^-»X) + t-c. 

+ ^{Hl*W^Hl- Hl*W^Hl) + ^. C . 

+ ^{QlB^d^ - iu*B^d„Ui + 2d*P' t d / X) + h.c. 

+ l -^(H:B^H u - H* d B»d^H d - L* 3 B»d„L 3 ) + h.c. 



ig Y 



{2e\B"d IJl e 1 + 2e* 2 B»d^i 2 + h.c.) . 



(55) 



C D + ^l 3 W+^P L ^l 3 ) + h.c. 



2 

ig Y 



(^B^Pl^+^B^P^IJ 
(-^ Hd B^P L ^ Hd * 2 Hd B^P L * 2 Hd ) 



- -f- (-^B^P**^ + 2^ dti B^P R ^> dii ) 
(2*? il B" 7 ^P fl *? il + 2*? i2 B" 7 „P fl *? i2 ) . 



(56) 



HereP L = (l+7 5 )/2 and P R = (1 -75) /2 are the left- and 
right-chiral projection operators respectively. \I/ U , are 
Dirac spinors representing the up- and down-type quarks 
respectively 



* — 



-ia 2 u* 



Qd,i 

-102(1*; 



(57) 



^] and tyf are Dirac spinors representing the neutrinos 
and charged leptons respectively (superscripts on Lj de- 
note the weak isospin component) 



L} 



Hi 



-ia 2 e* 



(58) 



and 



For a massless scalar 0, the rates for the decay of SU (2) w 
and U(l)y gauge fields to 4**<p pair are gj^M/967r and 
g Y y 2 M/96Tr respectively. Here M is the gauge fields mass 
and y is the hypercharge of <j>. 



The relevant terms for the decay of gauge fields to 











(59) 



For a massless fermion ip, the rates for the decay of 
SU(2)w and U(l)y gauge fields to -tpip pair are g^M/^&ir 
and g Y y 2 M/4:8ir respectively. Here M is the gauge field 
mass and y is the hypercharge of ip. 
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C. Fermion interactions 

The inflaton couplings to the fermions are found from 
the following part of the Lagrangian 

3 

C D V2gwYl \l{WIT{w 2 Li) + l\w\T \ia 2 L 2 ) 

i=l 

+V2g Y [elB t (ia 2 e 3 ) - iLjfl*(wT 2 ii) - ^S*(ia 2 £ 2 ) 
+h.c. (60) 

Here W\ , W 2 , W 3 and B are the gauginos of SU {2)w and 
U (l)y respectively. Superscript t denotes transposition, 
and 02 is the second Pauli matrix. 
The field content of L\, L 2l e 3 is 



are the supersymmetric partners of the W + and W 
gauge fields, see Eq. (Q2 



Li = 



L,= [ fj, e 3 = V 5 , (61) 



We find the following Dirac spinors from pairing -0's and 
gauginos (which are left-handed Weyl spinors) 



* 2 



* 1 




B 



(62) 



where 



W 



Wi - iW 2 
V2 



Wi + iW 2 
V2 



The ^F's decay into (s)quarks, the Higgs and Higgsino 
particles, and the L 3 , e±, e 2 (s)leptons. The relevant 
interaction terms for these decays are 



C d gw{QZ,i*iPL*d,i + Lj** x P L *? i3 + h.c.) 

+ g w (Hl*y 1 P L * 2 Hv +H 1 d *V 1 P L * 2 Hd +h.c.) 

+ 9w{Ql^2P L ^u,i + L 2 3 *^ 2 P L ^l 3 + h.c.) 

+ g w {Hl* V 2 PlK u +H 2 d **2P L K d +h.c.) 

+ ^(Q* u ^ 3 P L ^ u>i - Q* dti * 3 P L * d ,i + h.c.) 

+ ^(L$** a P L *l 3 - L 2 3 *W L *lz + h.c.) 

+ ^(H^ 3 P L ^ Hu -H 2 *<Sf 3 P L * 2 ^+h.c.) 

+ ^(H d *^ 3 P L ^ Hd - H 2 *^ 3 P L * 2 Hd + h.c.) 

+ ^r-(Q: 1 * 4 Pl*u, 1 + Q* d i^iPL^dA + h.c.) 


+ ^(1^^^3+1^*4^^,3+11.0.) 

+ ^y I (^u**4PL*k + Hi* 94P L 9% v + h.c.) 
X - //y (-H^Pl^ - H 2 *y 4 P L ^ 2 Hd + h.c.) 



2 

V2gv 
6 

V2g Y 



(-4Ui* 4 -PR*«,i + 2d i #4-F , R*<M + h.c.) 

(2ei* 4 PR*n + 2e 2 *4i > K*? 2 + h - c -) • ( 64 ) 



(63) 



The rate for the decay of ^1.2,3 to a massless scalar and 
its fermionic partner is g 2 v M/32ir, while that of ^4 is 
gyy 2 M /32ir . Here M denotes the mass of \& and y is the 
hypercharge of the scalar (and its fermionic partner). 
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